Projectile-shape dependence of impact craters in loose granular media 
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We report on the penetration of cylindrical projectiles dropped from rest into a dry, noncohesive 
granular medium. The cylinder length, diameter, density, and tip shape are all explicitly varied. 
For deep penetrations, as compared to the cylinder diameter, the data collapse onto a single scaling 
law that varies as the 1/3 power of the total drop distance, the 1/2 power of cylinder length, and 
the 1/6 power of cylinder diameter. For shallow penetrations, the projectile shape plays a crucial 
role with sharper objects penetrating deeper. 

PACS numbers; 45.70.-n, 45.70.Cc, 83.80.Fg 



Anyone who has walked across a beach realizes that 
even though the sand moves under the pressure of a step, 
at some point their foot is supported. Someone who is 
running leaves a deeper and larger footprint, or crater, 
than someone who is walking. This phenomenon applies 
to anything dropped into sand; no matter how energetic 
the impact, the sand will eventually stop the projectile. 
For an object of mass m that is dropped from rest and 
that moves downwards a total distance H (including mo- 
tion throughout the impact), the gravitational potential 
energy mgH is transfered from the projectile to the sand. 
Some of this energy is dissipated and some of it is con- 
verted back to gravitational potential energy associated 
with the size and shape of the crater. Fundamentally, the 
energy is expended by an average stopping force (F) ap- 
plied by the sand as the projectile moves a total distance 
d from initial impact to rest: 



{F)d = mgH. 



(1) 



Thus, a straightforward measurement oi d vs H pro- 
vides information about granular mechanics, a topic of 
widespread interest Q, • A full understanding of the 
stopping forces and the transfer of energy may enable 
a fundamental explanation for such impact-related phe- 
nomena as the complex sequence of crater morpholo- 
gies lE^IEI^ and the spectactular granular jets formed 
by a collapsing hole 0, II] ■ 

The penetration of spheres into loose granular media 
has been reported previously. In Ref. 9], the total drop 
distance H was varied by up to a factor of 1000; the 
projectile density pp was varied by a factor of 60; the 
projectile diameter Dp was varied by a factor of 4; the 
grain density pg was varied by a factor of 2; and the tan- 
gent of the draining repose angle p was varied by a factor 
of 2. The resulting penetration depths are reproduced in 
Fig. ^ along with the empirical scaling law 
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Evidently, the crater depth does not vary simply with 
either the kinetic energy of the sphere, ~ ppDp'^H, or 
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FIG. 1: Crater depth vs scaled total drop distance. All data 
are from Ref. for spherical projectiles. The solid line is the 
power-law scaling given by Eq. @. 



its momentum, ~ ppDp^ H^/^ . In Ref. Q a different def- 
inition of crater depth, suitable when the spheres com- 
pletely submerge, gave different scaling. 

The behavior in Fig. ^ cannot be explained by cur- 
rent theories. Granular hydrodynamics and a rate- 
independent "plowing" force both predict different scal- 
ing than Eq.|21(see discussion in Ref. Q)- But how gen- 
eral is the observation? In particular, how does the crater 
depth depend upon the shape of the projectile? One 
would expect a pointed projectile to penetrate deeper 
than a blunt projectile. Perhaps this may be accounted 
for by the value of the numerical coefficient in Eq. ||2Jl. 
Or perhaps the functional form of depth vs drop distance 
may be entirely different, with d ^ H^/^ arising as an ac- 
cident of the projectile's spherical shape. 

To explore the role of projectile shape on impact cra- 
tering, we perform a series of experiments with the same 
basic protocol as in Ref. Namely, we fill a 1000 mL 
beaker (diameter 4") with dry, noncohesive glass beads 
(diameter 0.2 mm, draining respose angle 24°). The 
depth of filling is 2" for most trials, but is increased to 4" 
for deep impacts. The density of the granular medium 
is Pg — 1.51 g/cc, which corresponds to random-close 
packing with a volume fraction of about 63%. As in 
Ref. 9], this state is achieved by swirling the sample 
horizontally, at first rapidly and then ever more slowly, 
until the surface is level and at rest. Instead of dropping 
spheres, we now drop a wide variety of cylinders and 
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FIG. 2: Crater depth vs drop distance for four different cylin- 
drical projectiles, all with diameter Dp = 0.25" and length 
L — 2". The solid lines represent fits to a 1/3 power- law, 



measure their depth of penetration. Cylinder materials 
include wood (poplar, pp — 0.5 g/cc), Aluminum (Al, 
Pp — 2.7 g/cc), and high-density Tunsten Carbide (WC, 
Pp — 17.1 g/cc). We use four different cylinder diameters. 
Dp e {0.25", 0.50", 1.0", 2.0"}, and four different cylin- 
der lengths, L e {0.5", 1", 2", 4"}. To vary the projectile 
shape, progressively sharper conical tips are machined 
at one end of these cylinders; we use four different cone 
angles, {180° (i.e. flat), 135°, 90°, 45°}. To achieve a ver- 
tical orientation at impact, the cylinders are suspended 
by a thin loop of thread, tied through an accurately ma- 
chined hole. For each cylinder, the drop distances are 
varied as widely as possible within certain limitations. 
The minimum is set either by the penetration depth for 
a cylinder placed infinitessimally above the surface of the 
medium, or by the tendency of (long or light) cylinders 
to fall over after a shallow penetration. The maximum 
drop distance is set either by the depth of the sand or 
by the tendency of cylinders to begin to tumble during 
free-fall. No trial is accepted unless the cylinder remains 
vertical after impact. 

Four sets of penetration depth vs drop distance data 
are displayed in Fig. |2 all for cylindrical projectiles of 
length 2" and diameter 0.25". This includes WC and 
Al rods with flat ends, as well as Al rods with 135° and 
90° degree tips. Evidently WC penetrates deeper than 
Al, in accord with their densities. The three Al rods all 
penetrate to approximately the same depth, in spite of 
their different tip shapes. Furthermore, all data are ade- 
quately described by a 1/3 power-law, d ^ just as 
for spherical projectiles. Therefore, the projectile shape 
seems to play no crucial role and the prior scaling law is 
not special to spheres. 

The results in Fig. |21 severely constrain the general- 
ization of Eq. ^ to aspherical projectiles. A natural 
hypothesis is that the crucial parameter is the mass per 
unit cross sectional area, {2/3)ppDp for spheres and ppL 
for cylinders. For cylinders, equal variation of density 
and length should lead to equal variation of penetration 
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FIG. 3: Crater depth vs scaled total drop distance. Data is for 
cylinderical objects, divided into graphs by cylinder diameter, 
as labelled. Values of L varried from 0.5", indicated by the 
smallest points, to 1.0", to 2.0", to 4" by the largest points. 
Tungsten Carbide is denoted by x, flat-ended Aluminum by 
□ , 135°-tip Aluminum by Q, 90°-tip Aluminum by V, 45°-tip 
Aluminum by and flat ended wood by ffl. The line in each 
subplot is Eq. @. 



depth. Thus, we must have 
d = 0.14 



1/3 



(3) 



for consistency with Eq. 

To test this expectation, we display depth data 
vs scaled drop height, [{ppL)/{pgp.'^)]^^'^Dp^^^H, in 
Figs. Ofa-d). This allows us to plot Eq. as a sin- 
gle curve, onto which we look for all data to collapse. 
Each subplot is for a different rod diameter; symbol types 
distinguish rods according to their length, density, and 
tip shape. Evidently, the degree of collapse is better for 
deeper penetrations and is worse for shallower penetra- 
tions. In Fig. for the 0.25" diameters, the quality of 
collapse is comparable to that found earlier for spheres. 
Comparison of the four subplots shows that the onset of 
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collapse appears correlated with rod diameter. In partic- 
ular, collapse by Eq. Q seems to be attained when the 
penetration is larger than the rod diameter. For shal- 
lower penetration, Eq. © fails and the actual depth de- 
pends on tip shape. In this regime, sharper tips cause 
deeper penetration. The tendency of the data to deviate 
from d ^ H^/^ for small H, concomitant with the onset of 
shape dependence, can also be seen in close re-inspection 
of Fig.El 

In conclusion, the crater depth vs drop height scal- 
ing observed previously for spheres is robust and can be 
simply generalized to projectiles of different shape. The 
power law is not an accident of projectile shape. 
The only accident is that it holds for spheres even when 
the penetration is not greater than the sphere diame- 
ter. The independence of deep penetrations on projectile 
shape, as observed here, should heighten the importance 
of finding a full theoretical understanding. It also offers 
a clue: The dissipation of the projectile's energy in the 
granular medium must occur throughout a volume that 
is much larger than the crater. Presumably, this involves 
a cooperative effect of many grains, such as long force 
chains successively created then broken by the projectile 
impact. 

To bring this manuscript full circle, we return to a 
walk on the beach. In particular, we generalize the depth 
scaling law to projectiles with non-circular cross section 
and we predict the depth of footprints in dry noncohesive 
sand {pg — 1.59 g/cc, p — tan(38°) Making three 

substitutions, ppL with the mass per cross sectional area 
m/A, Dp with 2^ A/n, and H with v'^ /{2g), we arrive at 



The area factors have been distributed to make it more 
apparent that this expresssion is dimensionally correct; 
at constant mass, the depth decreases with area as d ~ 
^-5/12 Pqj, inipact by a foot, we take A = 8 x 30 cm^, 
m ~ 70,000 g, and v ~ 100 cm/s (corresponding to a 
drop height of 5 cm). For these numbers Eq. predicts 
a footprint depth of 6.7 cm, in reasonable accord with 
experience. 
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